General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



(NASA-TH-8616 7) SiflOLTANEOOS OCEAN N85-16273 

CBOSS-SECUON AND BAINFALL EEASOBEHENTS PBOfl 
SPACE HITH A NADIB-PCINTING EADAH (NASA) 

46 p HC AOJ/MF A01 CSCL 04B Unclas 

G3/4 3 13002 

f\J/\SA 

Technical Memorandum 86167 


SIMULTANEOUS OCEAN CROSS- 
SECTION AND RAINFALL MEASURE- 
MENTS FROM SPACE WITH A 
- NADIR-POINTING RADAR 


Robert Meneghini and David Atlas 


NOVEMBER 1 984 


National Aeronautics and 
Space Administration 

Goddard Space Flight Center 
Greenbelt, Maryland 20771 


TM-S6167 


’...^ 


SIMULTANEOUS OCEAN CROSS-SECTION AND RAINFALL MEASUREMENTS 
FROM SPACE WITH A NADIR-POINTING RADAR 


Robert Meneghini and David Atlas 


November 1984 


GODDARD SPACE FLIGHT CENTER 
Greenbeit. Maryland 20770 


ABSTRACT 


In the case of a nadir^looking space-bome or aircraft radar in the presence of rain the 
return power corresponding to secondary surface reflections may provide information on the 
scattering properties of the surface and the precipitation. The object of the present study is 
to evaluate a method for determining simultaneously the rainfall rate and the normalized back- 
scattering cross section of the surface, o°. The method is based upon the mirror-reflected 
power, P^, which corresponds to the portion of the incident power scattered from the surface 
to the precipitation, intercepted by the precipitation, and again returned to the surface where 
it is scattered a final time back to the antenna. Two approximations are obtained for 
depending on whether the field of view at the surface is either much greater or much less than 
the height of the reflection layer. Since the dependence of P^^j on the backscattering cross 
section of the surface differs in the two cases, two algorithms are given by which the path- 
averaged rain rate and o° can be deduced. We also discuss the detectability of P^. the rela- 
tive strength of other contributions to the return power arriving simultaneous with P^^^, and 
the validity of the approximations used in deriving P^. 
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SIMULTANEOUS OCEAN CROSS-SECTION AND RAINFALL MEASUREMENTS 
FROM SPACE WITH A NADIR-POINTING RADAR 

1. INTRODUCTION 

The importance of measuring precipitation globally has been discussed on many occasions 
(e.g.. Atlas and Thiele, 1981). In the latter work, a broad spectrum of experts reviewed the 
state of the art and prospects for making precipitation measurements from space. Siifcsequently. 
Atlas et al. (1982) assessed a variety of spaceborne radar methods. This paper is a variation 
upon one of the methods discussed there, namely the surface reference method or SRT. The 
latter technique is treated in greater detail by Meneghini et al. (1983). One of its major short- 
comings was the need to estimate the normalized backscattering cross section of the surface. 
a°. independently in order that the rainfall rate could be measured. 

We shall show that for a nadir-looking radar, the present method offers the potential of 
measuring both the ocean radar cross-section a° and the path average rain rate simultaneously. 
This represents an improvement over the SRT method in that a° can be obtained directly 
from the measured powers. 

This paper is a revision of an earlier, preliminary study on the same subject (Atlas and 
Meneghini. 1983). Here we derive expressions for the mirror-rcllected power under fewer re- 
strictions and discuss the error sources in more detail. Nevertheless, a rigorous error analysis 
is beyond the scope of this paper. Therefore, the study should be viewed as a first step in 
addressing the questions of the nature of the secondary surface reflections and in determining 
whether meteorological information can be extracted from such data. 

2. BASIC CONCEPT 

Tlie basic concept of the method is an extension of the approach proposed by Eckerman 
et al., (1978). In essence, they proposed mea.surements of the surface radar cro.ss-section in 
the absence of rain from either measurements outside of the rain in the immediate vicinity, or 
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from a look-up table generated from prior observations. The path average rain rate would 
then be deduced from the microwave attenuation calculated us a difference between the echo 
expected in the absence of rain and that actually measured, and using well established relations 
between attenuation and rain rate. Dual wavelength versions of this method also have been des- 
cribed (Moore, 1981; Atlas et al„ 1982). 

The primary limitation of the method focused around the accuracy with which one could 
expect to estimate the surface radar cross-section. Indeed, the minimum detectable rain rate 
was found to be dependent upon the RMS error of this estimate. 

It is noteworthy that Inomata et al. (1981) conducted an aircraft experiment with two 
radars operating at 0.86 and 3.2 cm. By simply taking the ocean surface cross-section. 10 log 
a°, at nadir to be a constant 10 db at 0,86 cm. they were able to compute the attenuation 
at this wavelength and the average rain rate between the surface and the melting level. These 
values correlated well with those computed from measurements of the 3.2 cm reflectivity factor. 
Z, and resulted in a Z-R relationship close to those generally accepted (Battan, 1973), More 
recently, Fujita et al, (1984) have found fair to good correlations between path-averaged atten- 
uation as estimated from the surface reference method and a dual-wavelength algorithm (Fujita. 
1983). They note that the major source of error is probably caused fay the variability in 
Indeed the nadir ocean surface cross-section is sensitive to roughness and thus decreases with 
windspeed (Moore and Fung, 1979). This sensitivity has been exploited in estimating ocean 
surface wind speeds from space (Fedor and Brown. 1982). 

To reduce the errors caused by the variability in a° at least two approaches are possible. 
The first is to introduce a second radar wavelength and make use of the fact that the a° at 
the two wavelengths are well correlated, thus allowing measurements of differential path atten- 
uation. A second approach is to attempt to measure both a° and rainfall attenuation simultan- 
eously at a single wavelength. This paper is concerned exclusively with the latter approach for 


the special case of u nadir-looking radar over the sea. Although the method is applicable to sen- 
sing over the land as well, higher transmit powers are needed to obtain detectable signal levels. 

The scattering geometry is illustrated schematically in Fig. 1. Fig, la shows a nadir point- 
ing radar beam viewing the ocean through a storm. In the absence of rain, the echo power is 
a measure of the ocean backscatter cross section per unit area a°. With a rainstorm in the 
path, the power is attenuated by A db. which is a measure of the path average rain rate. R, 
(Atlas and Ulbrich, 1977). The latter authors show that, in the absence of vertical air motions, 
A(dB) is essentially linear with R regardless of the drop size distribution at 0.9 cm and approx- 
imately so in the band from about 0.7 to 2 cm. To get a true measure of A and thus of R 
(assuming the attenuating path length is known) it is necessary to know 

Figure lb illustrates the path of a ray from the radar through the precipitation, scattered 
from the ocean surface back to a slab of precipitation at height H. scattered back to the ocean, 
and finally reflected once more from the ocean surface and passing through the precipitation 
back to the radar receiver. The paths of the radar transmitter and receiver are displaced from 
one another for ease of displaying the incoming and outgoing rays. Of course similar mirror 
refiected echoes of the precipitation are received from all heights from which they are detec- 
table. 

Figure Ic shows a schematic of the radar return power as a function of range, or equiv- 
alently, time. The first contribution received is the direct backscatter from the precipitation, 

Pj, Ibllowed by the echo from the surface at height zero. P , At still later times the sum of 

o 

all scattering mechanisms corresponding to ray paths greater than (where is the satel- 
lite altitude) will arrive at the receiver. As indicated in the figure, one set of ray paths cor- 
responds to the mirror-refiected backscatter, Pjjj. In addition, other multi-path contributions 
will be present as well as the return from the surface received along the antenna side-lobes. If 
the mirror-retleetion component dominates, then it can be seen that a single series of measure- 
ments made along the range direction will yield P^, Pg and P^^. In the next section we .show 
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that under certain conditions the three unknown quantitites, a°, the rain attenuation and rain 
reflectivity can be found from the above measurables. 

3. DERIVATION 

To simplify the equations, we introduce the quantities Aj, j = 1, . , . . , n 

A: = k (s) ds ( i ) 

■ j o 

where Hj is the height of the jth range gate as measured from the surface. Letting the storm 
height correspond to the n**' range gate, i.e., = H„, then Aj, represents the one way 

attenuation through the precipitation. 

With this notation, the average echo power from the ocean surface under beam-limited 
conditions can be written (Appendix II) 

Pg = (Cg |o-0.2A„ (2) 

• where the symbols are defined in Appendix I, and the constant Cg is approximately 

Cg = 7 X 10-3 

Similarly, the echo power received directly from the precipitation at height Hj above the 
surface corresponding to a time of arrival of tj prior to that of the ground echo, is given by 

Pdf"V " ‘ 10-«'2 (A„ - Aj) 

where 

Cj = 3.5 X 10-3 P.J. L/0b“ (5) 

and where r?j is the rain retlectivity in the jth range gate above the surface. 

Under several simplifying assumptions, the expression for the mirror retlected power, 
measured at a time tj subsequent to that of the ground return can be written as (Appendix II) 
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( 6 ) 


P^(tj) = 3.5 X 10*3 


?J Vj L ff“ 10*°-2 (A„ + Aj) 


a” Pq 2 + 2.76 Hj2 
where is the radius of the field of view (FOV) given by 


Po 




(7) 


For notational convenience we let q = (a°l2J6 Thus, when the height of the 

reflection layer, Hj, is either much greater or much less than qp^, the mirror-refiected return 
power can be approximated by one of the following expressions 


C_' Tjj 0® r2 

,0-0.2 (A„ ^ Aj.) 


W = 




Cm 

JLJ to-0.2 (A„ ^ Aj) . u 


H. 


j « qpo 


( 8 ) 

(9) 


where F^ is the Fresnel reflectivity of the surface and 

= 1.26 X 10*3 \2 P.J, L 

Cn, = Cj = 3.5 X 10*3 X2 L / 0^^^ 


( 10 ) 

(ID 


The approximation to given by (8), apart from a constant factor that depends only 
on the particular antenna pattern function assumed, is the same as previously derived (Atlas 
and Meneghini, 1983). On the other hand, when Hj « qp^ and » Hj. Pj^ is the same 
as would be derived for the mirror reflected power from a perfectly smooth surface. In fact, 
in this latter case, a ratio of (9) to (4) yields P** 10**(-0.4 A) which implies that the ratio 
of the mirror-reflected power to the direct power near the surface, where Aj « 0, is approx- 
imately equal to the square of the Fresnel reflectivity of the surface. 

Tltere are several features of equations (8) and (9) that are of interest: (1) the precipi- 
tation parameters (reflectivity and attenuation) enter in the same way in both approximations; 
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(2) For Hj » qpg, is independent of the antenna parameters and satellite altitude but 
dependent on the height of the reflection layer, Hj, while for Hj « qp^, P,^ is inversely 
proportional to the area of the FOV but independent of Hj; (3) for Hj » qp^, P^^ is directly 
proportional to a® while for Hj « qpj,, P^ is independent of a°. 

Before discussing equations (8) and (9) it is worthwhile first to comment on why the form 
of Pjj, should change accoriding to whether Hj » qp^ or Hj « qp^. Consider a small volume 
element containing precipitation located along the vertical at a height Hj above tlie surface 
and upon which incident a ray scattered from the surface. If 2py » Hj then the majority 
of the energy scattered from the precipitation volume will be returned to the illuminated por- 
tion of the surface. However, if Hj » 2py, then of the energy scattered from the precipita- 
tion volume, only those rays nearly anti-parallel to the incident ray will be returned to the 
illuminated portion of the surface. In the former case, the precipitation is analogous to a mir- 
ror with a 'reflection coefficient’ of T7j', in other words, apart from the energy transmitted 
through the precipitation, all the energy incident on the precipitation is returned to the illumi- 
nated portion of the surface. On the basis of this difference we can e.\pect the form of Pj^ 
to change according to whether Hj is much greater of much less than the diameter of the 
FOV, i.e., 2pjj. Of course, this argument does not account for the fact that the actual quan- 
tity to be compared to Hj is not 2p^ but qp^ and therefore is a function of o° and as 
well as Py. 

An understanding of the dependence of P^ on a® may be gained by considering radiation 
incident on a surface area element which is scattered bistatically out to the precipitation, a 
portion of which is returned to a second surface element. The final process is the scattering 
from the secondary area element back to the antenna. The total mirror refiected power, there- 
fore. is given by the double sum of ail such primary and secondary scattering elements within 
the illuminated surface region. 
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In the notation of Appendix II the total mirror-reflected power is proportional to J where 

J = / / I(Xj, X2) e"^^ (1^1 1^ + 1^2^) dX| dX2 (12) 

and where X| , X2, are the locations of the primary and secondary surface scattering elements. 
The exponential term of the integrand represents the antenna pattern function at X) and x.,, 
with 0 = 7 = 0.86. It fnoj shown numerically that a good approximation to 

I(x,,X2)is 

l(Xj. X2) = Is " -2 (13) 

where, for a° » Hj » L 

Ij = jr L 0 “ / Hj2 (14) 

a-2 = ag2 = 8P2 Hj“ / o'* • (15) 

A heuristic argument for the form of Kxj. Xj) given by (13) is presented in Appendix II. 
In this approximation, I(X| . x,) is a function only of the separation between the scattering 
centers. As the separation distance goes to zero, 1(X], x>,) reduces to 1^ which is directly pro- 
portional to 0®. On the other hand, the distance between the scattering centers for which 

I (X[. X.2) = I5 e‘I 

is equal to a^ which in turn is inversely proportional to >/a® . We can interpret this last re- 
sult in the following way. about each primary scattering element on the surface. X|. there 
corresponds a secondary scattering region of dimension a^, centered at X| , which contributes 
the majority of the power returned to the antenna. This last statement must be qualified, 
however, if a^ » since in this case the radius of tile FOV, determines the effective 
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dimension of the secondary scattering region. In other words, the secondary scattering area 
can be defined approximately by 

Ag = min(jT a^^, 


To underst«in<4 this in a more quantitative manner we assume for simplicity that X|, 
centered at the surface area eiement AA(, is located at the center of the FOV, i.e. I X| | = 

Pl =0. The contribution to the mirror reflected power is then 

j| = AAj / Ij p dp 

0 


or 


Thus 


JT M, ij JJ^2 


J I “ 7T AA I 


: 7 » Po 

Ug ' 7 P() 


Since and are directly proportional to a®, the above equations show that J| is 
independent of a° for ya^ « p^ while J| is directly proportional to a° for ya^ » p^^. 

This accounts for the a° dependence in equations (8) and (9). The dependencies of P,„ on 
the antenna parameters and satellite and retlection layer altitudes can be found by summing 
the contributions from all primary surface scattering elements (Appendix II). 

4. ALGORITHMS 

Case I: Hj » qp^; Height of retlection layer much larger than diameter of FOV 

To obtain an algorithm that yields path attenuation and thus mean rainfall rate, 
we note that equations (4) and (8) are functions of pj thus allowing its elimination. 
Explicitly, 
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10-0.4 Aj 


( 16 ) 


^ Cm' °° (H. - H)>^ 

Pj(- (j) Cj Hj2 

To eliminate a°, (2) is used in (16). Taking the logarithm of the result yields 

2Aj-A„ = 5 |-Log[P^(tp/Pd(-tj)Pg] + Sj + S,} (17) 

where 

Sj = Log (H^2 (h^2 . h. 2) r2/Hj2] (18) 

S^, = Log(C„</Cd Cg) (19) 

When the range gate at the storm top is chosen, = H^, (17) reduces to an estimate 
for the total one-way attentuation. A„; 

A„ =■ 5 Log fP„(t„)/Pg(- t„) P,l + S„ + S„} CO) 

From (20), a measurement of Hj, and the k-R relationship, the path-averaged rain rate 
can be deduced (Atlas and Ulbrich; 1977). 

The estimate of attenuation given by (20) depends upon the calibration constants C^^', 

C(j. Cg through the term S^. If we assume that detectable levels of P^ can be obtained in 
light rain (where the attenuation is negligible), then the left hand side of (20) is approximately 
zero. This impiies that a best-fit value of can be obtained from measurements of and 
Pj near the storm top under conditions of light rain. Returning to (20) and comparing it 
with (2) shows that a° can be found from the equation 

= (Pg H„2/Cg) i0“-2 % (21) 

The vertical profiles of reHectivity factors, J?j, and the associated rain rates Rj can be ob- 
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tained from A„, equation (20), and measurements of Pjj(- tj) via the HitschfeJd-Bordan algor- 
ithm (see Lu and Hai, 1980; Meneghini et al., 1983). 

Case II: Hj « qpy; Height of reflection layer much less than diameter of FOV 

As in the first case, we begin by taking the ratio of to Pj where P^ is now given 
by (9): 


P„,(li) _ C„ r“ IH„ - Hj)^ 

Pd<- 'j> Cj Ho" 


Noting that then 



(23) 


Setting j = n in (23) again provides an estimate for the total path integrated attenuation, 
Aj,. Estimates of a° and the range-profiled rain rates are found by a procedure identical to 
that previously described. It is not necessary to assume that and Cj are equal, however. 

Since these quantities are range independent, it is dear from (22) that the interval attenuation 
between range gates k and j can be found. Solving for Aj|^ where Aj|^ = k(s) ds. r, > 
we obtain 






Pm(‘k> Pd(- tj) . 



(24) 


For Hg » the second term is negligible and equation (24) is essentially the same as 
would be obtained for a dual-wavelength algorithm (Ecdes and Mueller. 1971) if we interpret 
P,„ and Pj respectively as the attenuating and nonattenuating wavelength returns from the jth 
and kth range gates. In the dual-wavelength algorithm, the reflectivity factors are assumed to 
be equal at each range gate. i.e.. independent of wavelength. The analogous assumption here is 
that the reflectivity of the directly viewed rain volume is identical to that of the rain volume 
viewed via the first reflection from the surface (see section 6.3). 
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5. PRELIMINARY RESULTS 


T|je first objective is to determine how the mirror reflected precipitation echo behaves 
with respect to that directly scattered and that received from the surface. We also want to 
know the dynamic range of rain rates over which we can expect detectable mirror echoes, the 
smallest of the three measureables. For the purpose of these calculations we have assumed 
the radar» storm, and ocean param'eters shown in Table I. Initially we assume a very large 
antenna diameter of 7.S m to ensure that Hj » qp^ and therefore that equation (3) is ap- 
plicable. The effects of reducing the antenna size on are discussed below. 

Table I: Model Radar, Storm, and Ocean Parameters 


Parameter 


Value 

Wavelength - X 


1.87 cm 

Peak Power - Pj 


10 kw 

rulse Width — T 


1.33 sec 

Antenna Diameter - 

0 

7.5 m 

Radar Noise Figure 


7.0 dB 

Orbit Height - 


500 km 

Top of rain - H, 


3 km 

Ocean Cross-Section 

-0* 

10 dB 

Range Resolution - 

L/2 

0.2 km 


We aUo UK o* ■ 10 which corresponds to a wind speed of about 10 ms*^ Table 2 
values of 10 Log o* versus wind speed for the SEASAT 13.5 GHz radar (Fedor and 
«fown> mZi 


2= Normalized Backscatter Cross-Sections (10 Log o*) 
at "ifadir Versus Windspeed (Fedor and Brown. 1932) 

W^vpeed ims'’) i J 5 lO 20 50 lOO 

€>OSi«Se£l»flfdB) 105 13 5 118 9 7 7.2 5.1 4.0 

tkt feudts ^ in Fif. 2, AU vurv^i are signal to noise versus nun rate on a single 
|NAi« To 4 unifomt sale for all three quantities tlic plotted values of 10 Log(Pg/N). 


II 


corresponding to the sea echo, have been reduced by 60 dB while the values of the direct rain 
return to noise, 10 Log(Pj/N) have been reduced by 11 dB. Clearly, with such a powerful 
radar there is an abundance of surface signal with the signal to noise curve remaining approxi- 
mately constant at about 73 dB up to moderate rain rates beyond which the two-way attenua- 
tion through the 3 km depth of rain increases sharply. 

The curve of 10 Log(Pj/N), corresponding to the echo received directly from the precipi- • 
tation in the range gate immediately above the surface, exhibits adequate dynamic range with 
values exceeding 7 db over the full range from 0.1 to 100 mm/hr. At R < 10 mm/hr, 10 
Log Pj is proportional to 15.2 Log R corresponding to the assumed Z - R relation, Z = 424 
R*-^^. At R > 40 mm/hr, the curve is dominated by attenuation and the slope is proportional 
to “ 22 Log R corresponding to a two way path through rain with an attenuation coefficient k 
= 0.042 R*-* db/km. That is to say that a 10-foId increase in R produces a 22 db increase in 
attenuation and corresponding decrease in received power. 

For present purposes the most interesting curve is 10 Log (P^/N) corresponding to the 
range gate at the top of the rainfall at 3 km height. At low rain rates where attenuation is 
negligible P^^ is approximately 1 1 dB greater than P^ and the curves are parallel to one another. 
As the rain rate increases, the behavior of the ratio P^/Pj is dominated by the attenuation 
term for a two way path between the surface and height ’Hj. For the case shown, the gate 

for the mirror-reflected power is positioned at the top of the rain area, Hj = 3 km. so that 

the attenuation terms for P^^ and P^j are half that for Pj,,. When the range gates of the mirror 
reflection and the direct precipitation returns are near the ocean surface, the attenuation terms 
in Pjj^, Pj and Pg are approximately equal. 

As mentioned earlier, the size of the antenna was chosen to insure the validity of equation 

(8): i.e., using the approximations 0^ = \/2D and D = 7.5 m yield Py = 0.62 km 

which is much less than the maximum height of the reflection layer. For D = I m, « p^ 
so that P^ is given by equation (9). As shown in Fig. 3. P,„/N is reduced by about 7 dB from 
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the previous case and the range of rain rates for which P^/N exceeds unity is between about 
i.O mm hr'* and 30 mm hr'*. It is also worth noting that at light rains (before the effects of 
attenuation are evident) Pj exceeds P^ by about 4.4 db. This difference, which is present as 
long as Hj « qp^ results entirely from the fact that a ratio of equation (9) to equation (4) 
is equal to F'* I0**(- 0.4Aj), which for light rain rates and F^ = 0.6, is equal to - 4.4 dB. 

In Table 3 we present the range of rain rates for which P„,/N exceeds unity for selected 
values of antenna diameter and radar wavelength, for values of the peak transmit power of I 
kw and 10 kw. We see that the 1 kw system is inadequate with a 1 m diameter antenna. At 
3.2 cm, a combination of 1 kw and 3 m antenna diameter would serve well except at the 
light rain rates. On the other hand, with 10 kw, a 1 m antenna diameter provides adequate 
signal to noise over the lighter rain rates up to about 8 mm/hr at 0.86 cm. but quickly looses 
signal at the larger rain rates where attenuation is excessive. Thus, in order to achieve suffici- 
ent dynamic range in rain rate, we require a dual wavelength system. This is the case for any 
method, radar or radiometric, which relies on rainfall attenuation (Atlas et al.. 1984). We 
should note, nevertheless, that the condition required for the present method. P^^, » N. is 
significantly more stringent than the condition Pg » N, needed in the surface reference 
methods. 

Table 3: Range of Rain Rates (mm hr'*) for which P,„/N e.xceeds 

unity. The results in brackets correspond to Pj = I kw; 
otherwise Pj = 10 kw. 

Wavelength (cm) 




0.86 

1.87 

3.2 

Antenna 

) 

0.1 mm hf*‘+8 mm hr'' 

1 30 

4 ^ 70 

Diameter 

1 

(1.25 mm hr'*-*3.2 mm hr'* | 

[nonel 

[nonel 

(m) 


O.I-^IO 

0.3 -<* 40 

1 .3 -*• 100 


3 

1 0.25-^51 

1 1.5 2.51 

(6.0 -*• 64) 



0.l-*-l0 

0.2 40 

0.5 -»■ 100 


7.5 

(0.25 -<• 6.5| 

(1 25 1 

(2.5 -*• 80| 


6. SOURCES OF ERROR 

There are a number of difficulties associated with the method that must be resolved be- 
fore establishing its feasibility. Some of the difficulties are common to all meteorological ra- 
dar methods; e.g., sampling errors and errors caused by tluctuations in the drop size distribution. 
Other errors are, for the most part, characteristic of the method itself, Unfortunately, a rig- 
orous error analysis requires a knowledge of the bistatic cross sections of the ruin drops, a 
specification of the antenna pattern and the evaluation of a seven-fold integral. In the follow- 
ing sections, we discuss in a qualitative manner the error sources and indicate the circumstances 
under which they may be significant. 

6, } Multi-path and side-lobe contributions ^ 

There are a number of other scattering mechanisms that correspond to ray paths whose 
distances are equal to the mirror reflected path. Since these contributions arrive at the re- 
ceiver at the same time as P^, they constitute a source of error. Two kinds of paths are the 
dominant contributors to this error. The first is the direct backscattered power from the 
surface along the antenna side lobes. 

To find the direct side lobe contribution that arrives at the receiver simultaneously with 
P^, we set the path length of the mirror rellection equai to the path length of the side lobe 
contribution. This condition along with a knowledge of D/X and the antenna pattern function 
specify the side lobu(s) of interest. As an example, we assume that the primary rellector has 
a parabolic illumination function which produces an electric field in the far zone of the form 
(Silver. 1949) 

E = 8 (x)/x2 

where x = ;r D sin 0/\ and Jt is the Bessel function of the first kind of second order. For 
a mirror-refiection height of 3 km and for the parameters of Table 1 , the side lobe that arrives 
at the same time as is well beyond the 1 0th side lobe and therefore negligible in compar- 


ison to . If the antenna diameter is decreased to I m, the nearest side iobe of interest is 

the fifth. To obtain an upper bound on this contribution we take the maximum value of 
the gain within the fifth side lobe and the region of scattering to be equal to the urea of the 
annulus irHgL. This procedure gives a side-lobe to noise ratio, P^i^/N, approximately 83 dB 
down from the Pg/N curve shown in Fig. 3 which is well below the corresponding P„,/N curve 
for this case. In the algorithm for the interval attenuation (equation |24j), P^ is needed at 
two distinct range gates. When the height of the retlection layer, Hj, is near the surface the 
side-lobe contribution arriving simultaneously with Pj„ will be larger. For example, if Hj = 

0.5 km and D = 1 m then Pgi /N is down by about 35 dB from the plotted P„/N curve in 

o 

Fig. 3, As compared to Pn,(0.5 km)/N this is an unacceptably high level. To reduce P5L un 
additional 25 db the antenna diameter would need to be increased to approximately 1.5 m. 

The second contribution that arrives at the receiver simultaneous with P,^ can be under- 
stood by visualizing u ray path incident along the main beam (nadir direction) which scatters 
out of the vertical direction. This ray is then scattered by the precipitation and received along 

one of the antenna side lobes. The volume of interest is approximately that portion of the 

“shcir* of two confocal prolate spheroids which contains precipitation. The foci of the ellipsoids 

are the antenna and the center of the surface illuminated by the main beam (Fig. 4). 

To Investigate this multi-path contribution to the return power we begin by equating this 
path length to the mirror retlection path which gives 

where p is the radial distance measured from the axis of the main antenna beam (vertical) to 
the multi-path rain volume AV. We also define ^ as the angle between the ray directed from 
the surface to AV and the scattered ray from AV to the antenna (Fig. 4). Note that for 
» Hj. is nearly equal to the bistatic scattering angle 0 shown in the figure. 
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As belbrc Hj is the height of the mirror-reflection layer while z is the lieight of AV which 
ranges from ground level to the storm top. Table 4 gives values of p. ^ and n for three values 

of Hj: n is the side lobe nearest the main beam into which the scattered energy from.AV is 

received (an n =» 0 imples that the energy is received along the main beam of the antenna). 

Tile results have been tabulated for an altitude of = 500 km, where we have assumed \ = 

t.87 cm, z = Hj = 3 km and parabolic illumination of the retlector. 


Table 4 


= 500 km 

Hj(km) p(km) 


D = 1 m 

n 

D = 5 m 

D = 7.5 m 

0.5 2.7 

41,8“ 

0 

0 

1 

] 4 

53.3“ 

0 

I 

1 

3 8.4 

71.3“ 

0 

3 

6 


The large values of In Table 4 indicate that the rain volume of interest is outside the 
cone into which most of the incident power is bistatically scattered from the surface (see Fig. 
5). Moreover, since \p also represents the bistatic scattering angle of the rain drops, only a 
fraction of the energy incident on AV will be scattered towards the antenna. However, in 
four of the cases shown in the table, a portion of the multi-path contribution will be received 
along the main beam of the antenna (n = 0) which produces only a small reduction in the 
one-way antenna gain. Unfortunately, it is difficult to estimate the magnitude of the multi- 
path contribution without evaluating the bistastic cross section of the droplets and without 
accounting for the variations of ^ and p throughout the full range of z. Although crude 
estimates indicate that this contribution is smaller than P^, more detailed calculations arc 
needed. It should also be noted that a Doppler radar may provide a means of discriminating 
P,„ from side lobe and multi-path contributions (Atlas and Matejku, 1984). 


6.2 Bistatic Rain Reflectivity 

In the derivation of P,^ the bistatic rain reflectivity was replaced by the backscatter re- 
flectivity. The conditions under which the approximation is vuiid are not as stringent as might 
be thought since the pairs of surface elements that contribute most to the return power. , 

are separated by a distance less than the minimum of and 2p^ (Appendix II). Since is 

given by Hj P y/dlo° the range of bistutic scattering angles is given approximately by 


ff > ^ > T - 0 


where 


<tt = Tan*' 


3p2 Vt 

For example, if <f> = Tan*' |( — =-) ), then for o® = 10 dB, and P^ = 0.6 the minimum 

0 

bistatic scattering angle is approximately 180° - 35° = 145°. 

6.3 Mismatches Between Rain Scattering Volumes 

Because the ocean surface is roughened by the wind and waves, energy incident on the 
surface is scattered in ait directions. This is why the mirror precipitation echo must be inte- 
grated over all scattering angles. The bistatic scatter cross-scction corresponding to i£q. (A9) 
after Ishimaru (1978) is shown in Fig. 5 with windspeed as a parameter. In order to assign 
a wind speed to each curve, we matched the value of 10 Log o° for direct backscatter (i.e., 

0 deg) to the corresponding values read from the nadir values of o° versus wincl.speed given 
by Fedor and Brown (1982) and reproduced in Table 2, Since the latter data were based on 
observations taken at vertical incidence, while the curves are theoretical, one should not regard 
them as rigorous. Nevertheless they indicate reasonable behavior. In essence they show that 
as the wind Increases the nadir surface cross-section decreases and that increasing energy is 
scattered to angles father off nadir. The inset table in Fig. 5 lists the full angle around nadir 
for which the cumulative power wliicli is scattered from the rain is equal to half the total 
power radiated back toward the radar. For example, at a windspeed of 10 ms*' (correspond- 
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ing to 10 Log a° 10 dB), we must include the energy scattered from withu' a cone of 34.4 
deg to account for half the total mirror reflected power at the iradur. 

As indicated in Fig. S, even light winds will roughen the ocean sufficieiuly to spread the 
retlected beams over fairly sizable angles. This has the effect of increasing the horiz'^^.ntal 
dimensions of the volume of precipitation which is viewed in mirror fetlection rciative to that 
which is viewed directly by the incident beam, the more so the higher the votuine element. 
The effects of beam divergence cail into question two assumptions made previously. In the 
derivation of it was assumed that most of the scattering takes place along the vertical, 
implying that the general path attenuation could be replaced by the attenuation along the 
vertical. However, as the beam divergence increases, the attenuation must be considered as a 
weighted average along the various directions of the surface scattered power. The errors re- 
sulting from the vertical scattering assumption can be studied quantitatively by including the 
general path attenuation into the definition of the integral l(X|, X 2 ) (equation (AI3)). A 
study of this type, however, is beyond the scope of this paper. Secondly, in the algorithms 
for path attenuation, we Implicitly assumed that the precipitation within the pulse volume 
seen in mirror retlection be identical to or of the same rellectivity to that viewed directly. 

Tliis can be a demanding requirement at a windspeed of 10 ms**, with the half angle of 17.2 
deg, and a pulse volume at 3 km height. The horizontal dimension of the mirror scattering 
volume is then increased by approximately 0.93 km on either side of the incident beam. For 
D = 7,5 m, \ = 1.87, = 500 km the diameter of the incident beam near the surface is 

approximately 1,2 km. In this case the mirror reflected scattering volume, at a height 3 km 
above the surface is over a factor of 6 greater than the direct rain scattering volume. On the 
other hand, for an antenna diameter of I m the ratio of the mirror reflected volume to the 
direct volume decreases to about 1,44. 

7. SUMMARY 

The method proposed here for the remote measurement of rainfall rate either from space 
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or uircrul't is motivated by the need to circumvent one of the limitations of the single wave- 
length surface reference technique; namely, that estimates of a” are required to extract meteoro- 
logical information from the surface return power, To provide additionai information on tiie 
surface and meteorological parameters, the mirror-refiected power, P^, is considered us a pos- 
sible candidate. corresponds to the portion of the power scattered by the surface, inter- 
cepted by the precipitation and returned to the surface where it is scattered a final time back 
to the antenna. Under a number of simplifications a closed form expression for P^ is ob- 
tained. The dependence of P^ on the backscattering cross section, a”, was found to depend 
upon the magnitude of the quantity Hj rA/a“ p^, where is the Fresnel reflectivity of the 
surface, Hj is the height of the reflection layer and.p^ is the radius of the FOV. Wlien tills 
quantity is much greater than one. P,^ is directly proportional to a°*. if it is much less than 
one, P^ Is independent of o'*. Algorithms based on these limiting cases are proposed for the 
estimation of path attenuation, rainfall rate and a°. 

Tile levels of P^ to noise power at nadir over the ocean indicate that detectable signal 
levels can be obtained from satellite altitudes over a fairly wide rain rate interval. This remains 
true even with a modest antenna dimension of I m. Thus, the method may provide an im- 
portant auxilliury measurement scheme for a space-borne meterologicul radar for tlie estimation 
of ruin and the buckscattering cross section of the surface. Estimates of this latter quantity 
will in turn permit one to investigate the dependence of the mean square surface slope on 
wind speed in the presence of rain. When tlie peak transmitted power is reduced from 10 kw 
to 1 kw, however, the range of detectable rain rates is reduced considerably. Furthermore, 
over land the decrease in P^ and o® relative to the ocean background indicates tliat even at 
10 kw. the mirror reflected power us measured from a low earth orbiter will be smaller than 
the receiver noise power for most rain rates and radar wavelengths, 

In the derivation of Pj„ a number of simplifying assumptions were made both to make 
the problem tractable and to understand the approximate dependence of P|„ on metcrological 


f 
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und surface scattering parameters. Although a qualitative discussion of the error sources is 
given, quantitative results are needed to check the . range of validity of the expressions for 
and to compare the relative magnitudes of side-lobe and other multi-path contributions that 
arrive at the antenna simultaneous with P^. 


FIGURE CAPTIONS 


Fig. I; Schemutic of the mirror-reflected echo. The mirror-reflection echos which arc 

returned subsequent to the surface return In Ic occur by means of a double reficc-. 
tion from the surface as shown in Ib. 

Fig. 2 j Curves of the signal to noise ratios for returns from the ocean, Pg. the precipitation 
viewed directly at a range Just above the surfaces Pj, and the mirror-rellected re- 
turn, P^, from a 3 km height. Radar parameters are given in Table 1. 

Fig. 3: Same as Fig, 2 except for an antenna diameter of I m. 

Fig. 4; Schematic of a scattering path that arrives at the receiver simultaneous with the 

mirror-reflected return. 

Fig, 5: Bistatic cross section of the ocean for nadir incidence. Inset table shows the conical 

angle from which 50% of the total power is received. 

Fig. All Scattering geometry of a ray incident from the radar to the ocean surface element 
M|, scattered to the precipitation volume AV, and returned to the surface element 
where it is scattered a final time back to the radar. 

Fig. A2; Schematic of the precipitation volume (ellipsoidui shell) that Contributes to the * 
mirror-reflected power for two fixed surface scattering elements. 
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APPENDIX 1 


TaWe of Symbols 


X 

- 

Wavelength 

h 

- 

Peak transmitted power 

Pl 

- 

Transmitted power waveform 

t 

- 

Transmitted pulse duration 

L/2 

- 

Range resolution 

Ho 

- 

Orbit altitude 

«i 

mm 

Height of rain layer 

H, 

- 

Height of storm 


- 

Ocean crosS'Section at nadir 


- 

Ocean normalized bistatic cross*scetion: incident wave in the T 



direction, scattered wave In the direction u. 

h 

- 

3 db beamwidth measured from the peak to the half power p< 


~ 

Effective antenna area 

% 

“ 

Antenna gain along beam axis 

GjiTi 

tat 

Antenna gain pattern along the ith direction on tninsmic 

g,iT> 


Antenna gain pattern along the ith direction on receive 



Bistalic rain reflectivity 

n 


Baekseattered rain retiectiviiy 

a 

tm 

Palh*averaged rain rate 

k 


«\ttenuaiton coefneient 

z 

*=Sr 

a 

Ram rcllectivity factor 

A 

«» 

Paih'ifltegrated attenuation 


sa 

Xlirror reneeted echo power 



Table of Symbols 

Pjj - Echo power directly scattered from precipitation. 

P_ - Echo power directly scattered from ocean surface. • ij 

6 * P 

N - Receiver noise power | 

c - Speed of light in free space i] 

IJ 



APPENDIX 11 


The purpose of the appendix is to derive approximate expressions for the return power 
P^. corresponds to the portion of the power scattered by the surface which Is intercepted 
by the precipitation, returned again to the surface and Is scattered back to the antenna. A 
general ray path contributing to P^ is shown in Fig. Ai and can be written as 




■*.(«- 1 n /'=> 

(47t)‘* (rg r, rj r 3 )“ 




* (abCtjttj) AA.) • A, (Al) 

where Q is that portion of the ray path within the precipitation and where the .symbols in {AD 
are defined in Appendix 1. 

Assuming that AA|. AAt, and AV can be converted into differential area and volume 
elements, and using the relation A^ = Gj^ir then P^ can be written 




(4jt)' 


I 


c,(fo) 0,(^3) 

2 2 
*•0 ^3 


dX| dx. 


X 






io-«-i dvj 


(A2) 


where the volume integral is taken over the precipitation, and dX|. dx, are elements of sur- 
face area. The range of integration on both X[, x, is the portion of the surface illuminated 
by the main beam. Let the transmitted pulse be rectangular and of duration r with P-p the 
peak transmitted power 


P, (t) = 


P.p [U(t) - U(t - r)l 


(A3) 


i 


where U is the unit step function. Therefore, 

3 

P,(t - S Ti/c) “ Pj (U(t - T - (r, + r2)/c) - U(t - T - (r, + r,)/c - r)] {A4) 

If Hg » pQ where is the spacecraft altitude and p^ is the radius of the footprint, 

then 

T =« 2Hq/c + (|x,l2 + Ix,^) /2H„ c (AS) 


where |X|I, IX'^I are the distances measured from the center of the footprint to Xj , x.,, respec- 
tively. Thus Pj is non zero for 

2H < r, + r, < 2(H + L/2) IA6) 

where L/2 is the ramie resolution, equal to cr/2, and 

2H = ct - 2H„ - (Ix,|2 + ix2p)/2H^ -CT (A7) 


The condition given by (A6) specifies the region of precipitation that contributes to P,^. 
In particular, for a fixed observation time and fixed Xp‘ x.^ the volume integral taken over 
the precipitation is that region contained within a prolate spheroidal ‘shell’ with foci at X| . x., 
(Fig. A2). 

Using a coordinate system whose origin is mid-way between X| , x., and whose x-axis 
passes through both points, the family of prolate spheroids within this region can be expressed 
by the equation (Morse and Feshbach. 1953) 

n 


f,nin ^ f ^ fmax 


(y- + z^) 

+ L = I ; f > |3 > a 

. 


(A8) 


^min 

^max 


(H“ + 

[(H + L/2)2 + or]'^' 


(A9) 
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The distances between the points X| and ^2 Siven by 2 Xq where 


Xg = 


To simplify the expression for P^, we make the following approximations 

(a) The bistatic rain reflectivity factor, t7(7'| *, 72 ), can be replaced by the backscattcred 
rain reflectivity, rj. 

(b) The majority of the power backscattered from the surface is directed aiong the ver- 
ticai direction so that the attenuation can be approximated by the attenuation along 
the vertical direction, i.e. 

/ k ds 2 /^s k dz + 2 k dz 
Q 0 0 

where Hj is the storm height. We assume, moreover that H can be interpreted as the height 
of the reflection layer and that H is nearly independent of X|. X 2 - 

(c) The spacecraft or aircraft altitude is taken to be much greater than the radius of the 
main-beam footprint; Hg » Pg. This implies tliat the ‘ray’ paths Xg, r^ siiown in 
Fig, (Al) are nearly vertical so that 


7,) = ff5|) 


where the angles 9 ^, 9 -, are defined in Fig, (A2). 

(d) The bistatic cross section of the surface for incident energy along the vertical which 
is scattered at an angle 9 is given by the scalar Kirchhoff approximation for a very 
rough surface (Ishimarir. 1978); 


f(0) 


2mf^ -•^tan“fi/2 

e ^ 


(I + cos 


(AlO) 


where is the Frcsnei reflectivity of the surface and m = where £ 

is the correlation length of the surface and |3 q Is the rms surface height. 

(e) Assuming a Gaussian antenna pattern, the antenna gains are written 

Gj(fo) = Gy e*<7^l/®b)^ » Gy 

Gflfj) = Gy Gy q-(^^P2' 


where is the beamwidth measured from the maximum to the half power point and where 

Gy = ff2/(20b)2 

7 = 0.831 

P =• y/O^H^ ^ 7/Po 


Tlie quantities pj , P2 respective distances from the center of the Held of view 

(FOV) to the scattering centers X| , X 2 on the surface (i.e. pj = |X||, p-, = Ix^l) while Py 
is the radius of the FOV. 

The validity of these approximations Is discussed in the text. With these assumptions 
(A2) can be written 




Gy2 Pj TJ 

(4jt)^ Hy-* 


•0.2(A 
10 " 


+ A ) 


(All) 


where 


J = / / 

+ Pn') 

l(2i| . Hi) tl^i 

(A12) 







f(0i) ffO-,) 



J = / 

1 dV 

(AI3) 


V 

r,- r 2 - 



A„ = k dz-. 

A = k dz 

(AI4) 


0 

u 


where the observation time t 

corresponds to the 

rellection layer height. 

H and where f(0) is 


given by (AlO). 
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To express the volume integration of (A 1 3) as an iterated integral, we first note from 
Fig. A2 that cos B-y = z/r-?, cos 6^ = z/r| and therefore cos ~ fi ‘-’os 0^/r-y. 

For any point (x, y, z) within the spheroidal ‘shell’ shown in Fig. A2 we have 

r, = ((x + X(j)2 + y2 + z2|'/j 

= ((X - Xjj)2 + y2 + z2)'/: 

where as before the origin of the coordinate system is taken midway between the two scat- 
tering centers where the x-axis passes through X| and Xj- Titerefore. 

Ti = (r,* + 4 xXq)*^^ (A 1 5) 

Choosing a spherical coordinate system with an origin at x = x„, an arbitrary point (x, y. 
z) in the reference systetn can be expressed by 

X = Xq + rT sin Ot cos 

y = Tt sin sin 0^ (A 16) 

z = r-» cos 0-y 

Substituting (A 16) into (A15) gives 


where 


so that 


r, “ T 2 F 


c r , ‘"‘o • « 1 

F = 1 + sin 0 -) cos + j 

L '2 “ " '-2“ J 


cos = cos d-)/F 


(AI7) 


(AI8) 


(A19) 
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Using (A 1 9) and (A 10) the function f(0|) can be written in terms of to. 
plicitly, f((?|) a f(r 2 , 02 ' ^ 2 ^ 


t(r->t 0')* ~ 


(I + cos 
Tlierefore {A 13) becomes 


exp ( - - (t - cos- 0i/F)/(I - cos 0,/F)2 1 (A20) 


Kxi.x-.) = 


sin 0 dr d0 dd 


(A21) 


^■0 0"O r«r^ r- F“ 

The limits on r cun be found by substituting equations (AlO) into (A8) and solving T 2 
when f = f^jn, Procedure gives 

rA(0» “ Xq^)/(H + Xg sin 6 cos 0) 

fB(0, 6) = ((H + L/2)2 - Xo^l/KH + L/2) + x^ sin 0 cos 0| (A22) 

In spite of the simplifications made, the integral (A12) with I given by (A2I) is formid- 
able, To further simplify (A2I) we use the fact that a good approximation to l(X|, X 2 ) is 

Kip i2^ ■ *s ~ ^2^“ 


where 


«s = 


7T L 

H 


0 ? = 1/a 2 = a°/8 h2 


IA23) 

IA24) 

(A2S) 


and where a° is the normalized buckscattering cross section of the surface for nadir incidence. 

To understand qualitatively why Kxj, x, ) can be approximated by this function we first 
note that us xt “*■ ip (A21) becomes 


ir/2 H+L/2 

Ki] < i] ) = Ij ~ ^ ~ — ^ 

0=0 r=H 


(A26) 
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Making the change of variable u = tan^(<?/2) and using (A 10), (A26) can be integrated 
to give 




I, = - (H + L/2)'' ) 


m 




-e-'")) 

li?—] 


(A27) 


From Ishamaru (1978), we have 


m * 2o®/r^ 


Since for most sea states m » I, we can negiect terms of order m'^, m*^. and e'*". Further- 
more, assuming that H » L/2 then Ij reduces to (A24). 

The second feature of l(Xj , >< 2 ), i,e„ its dependence on |X| - x,,| and a^ can be under- 
stood by noting that the integrand of i(X|, X2) is proportional to f(0| ) f(d,) where 

ffdj) a exp [ - -^(tan^ (0,/2) + tan’ (OnlZ) )| 

Letting 0j = ^ 2 = 0 we obtain f"(6) « e“* when tan 6 = 2>/m7(ni - 1), The angle 6 
cun be related to the distance between the scattering centers, a^. by noting that 


ag = 2H tan 9 

or, using the above result for tan 9 and m = 2a®/P^ » I, then 

ag = 4HF / 

Returning to (A23) and substituting this equation into (A 12) and expressing X| , k-, in 
cylindrical coordinates we obtain 


J = I, r r e-<“‘ + + P2-) 

P|=0 P2=0 




Pi 


d0| d^2 [ dp, dp2 


f y’’’ ^J20!“P|P2 cos((^| - 


(0| = -7T -7T 


(A28) 
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w 


To perform the integrals over ninke the transformation 


u => - ^2 

V = 0, 


The double integral over , 02 then becomes 


I, = / / ^ 1^2 cos u 

V«-T U«'JT-V 


^pcosu du » 2 / ePCosu dit » 2ff lo(p) 
u«-v-jr 0 

where 1^ is the modified Bessel function of the first kind, zeroth order. 
Therefore, 


•l = P|P2> 


Substituting <A29) into (A28) gives 


J = 4it^ *s / / PjPt) ^ * ^2 dp|dpo (A30) 

p2=0p,=0 


Carrying out the remaining integrations (Gradsiiteyn and Ryzhik. 1965) yields 


[2«2 + (J2j 


(A3I) 


Substituting (A3I) into (Al I) and using tiie definitions of a, ;3 and gives equation (6) 


of the tc-Kt. 


To eulculatc the ground return we assume that beam-limited conditions hold and we 


choose, as before, a Gaussian antenna gain. We obtain 


Pj Gy2 a» \2 


p = I0-O.2 A„ 

* 2’ »= (J2 H„" 
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Using Gy = definition of ( 3 , yields equation (2) of the text. The expres- 

sion for the direct rain return, equation (3) is essentially equivalent to that derived by Probert- 
Jones (1962) for a rotationally symmetric antenna pattern. 
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Figure I . Schematic of the mirror-renected echo. Tlie mirror- 
retlccted echos which are returned subsequent to the 
surface return in Ic occur by means of a double rc- 
(lection from the surface as sliown in !b. 
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Figure 3. Same as Fig. 2 excepj for an antenna diameter of I 










Figure 5. Bistatic cross section of the ocean for nadir 

incidence. Inset table shows the conical angle 
from which 50% of the total power is received. 


Figure Al. Scattering geometry of a ray incident from the radar 
to the ocean surface element AA| . scattered to the 
precipitation volume AV. and returned to the surface 
element AA-) where it is scattered a final time back 
to the radar 
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